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Abstract. Electron-impact Vi — > Vf vibrational excitations cross sections, involving 
rovibrationally excited N 2 (wi, J) and NO(vi, J) molecules (fixed J), are calculated for 
collisions occurring through the nitrogen resonant electronic state N^"(X 2 n g ), and 
the three resonant states of nitric oxide NO _ ( 3 S - , 1 A, 1 Complete sets of cross 

sections have been obtained for all possible transitions involving 68 vibrational levels 
of N 2 (X and 55 levels of NO(X 2 I1), for the incident electron energy between 0.1 

and 10 eV. In order to study the rotational motion in the resonant processes, cross 
sections have been also computed for rotationally elastic transitions characterized by 
the rotational quantum number J running from through 150. The calculations 
are performed within the framework of the local complex potential model, by using 
potentials energies and widths optimized in order to reproduce the experimental 
cross sections available in literature. Rate coefficients are calculated for all the 
(t>i, J) — ¥ (vf,J) transitions by assuming a Maxwellian electron energy distribution 
function in the temperature range from 0.1 eV to 100 eV. 

All the produced numerical data can be accessed at 
http://users.ba.cnr.it/imip/cscpal38/phys4cntry/databasc.html. 
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1. Introduction 

Vibrational transitions, induced by electron impact in diatomic molecules, play a key 
role in affecting the chemistry and the physics of many molecular plasma [lj. The most 
efficient low-energy process, leading to molecular vibrational excitation/de-excitation, 
occurs through a resonant mechanism, according to which the incident electron is 
momentarily captured by a neutral molecule, in general vibrationally excited, with 
the subsequent formation of a resonant negative molecular ion. Within a given time 
interval, the electron can acquire all, or part, of its original kinetic energy and may 
escape from the molecule giving rise, correspondingly, to elastic or inelastic vibrational 
excitations [2]. 

We will focus our discussion on electron-impact vibrational-excitation collisions 
involving N2 and NO molecules. The interest in these species relies on their relevance, 
as plasma components, in several fields of scientific and technological importance: from 
planetary atmosphere, including environmental and spacecraft re-entry problems, to 
nuclear fusion and industrial plasmas [31 Hj . 

Resonant vibrational excitations (RVE), for these two molecules, occur via to the 
following mechanisms: 



In both cases, the two molecules are initially in their respective ground electronic state 
and in one of the associated vibrational levels, labeled by the quantum number Vi, 
ranging from through 67 for N 2 and 54 for NO and a fixed rotational quantum number 
J. The capture of the incident electron leads to the formation of a resonant negative ion 
which, after autoionization, decays back to some level v / of the neutral molecule. NO - 
contributes to process ^ with the three resonant electronic states 3 £~, 1 A, 1 S + , while 
process ([!]) proceeds through the resonant X 2 U g state, responsible for the well-known 
2.3 eV resonance peak in electron-^ collision experiments [5]. 

Experimental measurements 0, [7J [8j [TUl HTj [12] and theoretical calculations [13, 
dU [TSJ [TBI CEO HS1 [13 ED] of cross sections exist in literature for the above processes. 
However, despite the broad production of cross sections data, especially for nitrogen, 
all the previous works are restricted to excitations between the first few vibrational 
levels. This circumstance is particularly limiting for models of those plasma systems in 
which strong non-equilibrium conditions exist. In these plasma systems the vibrational 
distribution is no longer Boltzmann, so that the high vibrational levels can become 
over-populated. An appropriate example of non-equilibrium conditions is provided by 
the air plasma generated by the impact of a supersonic space vehicles during its re-entry 
phase in the earth atmosphere. During the air-vehicle interaction, the kinetic energy 
transferred to atmospheric components induces complex physical and chemical processes 
and may, in particular, activate the internal degrees of freedom of the molecular species 
whose densities can strongly deviate from Boltzmann distribution. A realistic model, 
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aimed to predict the chemical and physical evolution of these plasmas, relies therefore 
on a detailed description of the state-to-state processes occurring at the molecular level, 
with particular consideration of the electron-molecule energy exchanges, whose role is 
of fundamental importance in determining the level population of the chemical species. 

In this paper we present complete sets of cross section data for both the excitation 
processes ([TJ and ([2]), connecting all the vibrational levels of N2 and NO molecules, 
and for an incident electron energies interval below 10 eV. We have also extended the 
calculations to rotationally elastic resonant excitations of both the molecules (J < 150 
and A J = 0). The cross sections are evaluated by resorting to the local complex potential 
(LCP) model as formulated in Ref. [2T]. The LCP approximation is derived from the 
more rigorous nonlocal model [22J, which is presently the best available theory for the 
description of the resonant collisions. Comparison of the two models, however, reported 
for N2 and NO in Refs. [15] and [H], respectively, for some transitions, shows that the 
results of the two methods are in good agreement, and the observed differences, although 
significant for testing the models, are quite inconsequential for any practical purpose. 
This is a favorable circumstance from a computational point of view, as the LCP model 
is less demanding, in terms of computer times, compared to nonlocal counterpart, 
especially if we consider, for example, that the complete sets of vibrational data include 
68x68 and 55x55 possible transition cross sections to be evaluated for each energy 
grid point. Beside the non negligible time-saving features, the LCP approximation also 
provides quite accurate results demonstrated, as shown later, by the good agreement of 
the calculated cross sections with the experimental measurements. 

Cross sections provide the basic input data for a plasma kinetic model, from which 
the rate coefficients can be obtained once the electron energy distribution function, for 
the particular conditions of the plasma system, is known. In many situations, however, 
when the electrons can be considered thermalized at a given temperature T e , the rate 
coefficients can be obtained by convoluting the cross sections with a Maxwellian energy 
distribution function, according to the equation: 

I g / 1 \ 3 / 2 r°° e 
K ^- (T ' )= i^rAwJ l/EEe-^^iE), (3) 

where m e is the electron mass, k is the Boltzmann constant, e#, = Vf — Vi is the threshold 
energy and a Vi ^ Vf (E) is the cross section for the Vi — > Vf resonant vibrational excitation 
for electron energy E. Using the calculated cross sections for both processes ([TJ and 
fl2J), we have also calculated the corresponding rate coefficients, by Eq. (|3]), for all the 
transitions and in the range of electron temperature from 0.1 eV up to 100 eV. 

The organization of the paper is as follows. In the next Section we briefly 
describe the main equations of the LCP model and give some computational details. 
In Section [3] we compare the calculated cross sections with some available theoretical 
and experimental results of the literature, and in Section [4] we illustrate few particular 
examples of the calculated cross sections and rate coefficients. In Section |5j finally, some 
conclusive comments will be provided. 
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2. Theoretical model and numerical details 

In this section only a brief account is given of the relevant equations of the LCP model. 
A comprehensive formulation of the resonant collisions can be found in [21J and in 
references therein. 

The cross section for a given rovibrational (fj, Jj) — » (vf,Jf) transition, occurring 
through an electron- molecule resonant collision is given by |21j : 

. , Q r 2 7T^ TV? kf I ,2 , .v 

= ^-j^ Kx/i^lOl , (4) 

where p r and (7 represent, respectively, the product of the degeneracy factors and 
spin-multiplicity of the resonant and neutral electronic states. The equation includes 
a multiplicative factor of 1/2 in order to account for the spin- multiplicity of the 
unpolarized incident electron. The initial and final electron momentums are defined by 
kf = 2m e €i/h 2 and k 2 = 2m e €f/h 2 , where e« and e/ are the initial and the final kinetic 
energies of the free electron. x v .j (R), is the wave function of the final rovibrational 
state of the target molecule and £(R) the nuclear wave function of the resonant state, 
both depending on the internuclear distance R. The exit amplitude (R) is the value 
of the interaction matrix element, Vk(R), coupling the discrete and continuum states of 
the colliding system [21] for the exit channel k = kf. 

Once the coupling matrix Vk(R) and the potential energy curve of the resonant 
state, V~(R), are known, the wave function £(R) can be obtained in the LCP model by 
solving the nuclear resonant equation: 



h 2 d 2 JiiJi + ^h 2 , T _ , D , , . , D , i 



+ o m + V ~( R ) + A ( R ) - ~ T ( R ) ~ E 



2fidR 2 2/iR 2 v ' v ' 2 

= -V^R) Xv^iR) , (5) 

where \i is the nuclei reduced mass and E = e« + E Vit j. — e/ + E Vft j. is the total energy. 
E v j. are the rovibrational eigenvalues of the target Schrodinger equation: 



k 2 d 2 J t (J t + l)h 2 



Xv ,j t (R) = E v ,jXv,jAR) ■ (6) 



2fidR 2 2/ii? 2 

Here Vq(R) is the electronic potential energy function of the neutral molecule. In Eq. (J5|, 
T(R) and A(R) are respectively the resonance width and the level shift, and finally, 
Vk^R) is the entry amplitude describing the capture of the incident electron by the 
molecule. 

The potential curves, V~(R) and Vo(R), needed in Eq. (JsJ) and ^ respectively, 
have been expressed Morse-like function: 



U(R) = D e [l-e- a(R - R ^] 2 + W , (7) 

where D e is the dissociation energy for the considered electronic state and R e is its 
equilibrium bond-length. W is the zero-point energy and a is a constant. The Morse 
potential parameters for the nitrogen molecule were obtained, by a fitting procedure, 
from the Gilmore potential [23] • The dissociation energy, measured from the potential 



Resonant VE electron-^ and -NO cross sections 



5 



Table 1. Morse parameters for N2, N 2 , NO and NO potential states (Eq. ([7])) and 



reduced mass /1 for the neutral species. 





N 2 (X x £+) 


N 2 -(X 2 U 9 ) 


no(x 2 n) 


NCT( 3 £-) 


NO"( 1 A) 


NO"( 1 E+) 


D e (eV) 


9.906 


8.638 


6.61 


5.161 


5.411 


4.95 


R e (a.u.) 


2.069 


2.230 


2.175 


2.393 


2.38 


2.37 


a (a.u.) 


1.428 


1.20 


1.48 


1.20 


1.18 


1.20 


W (eV) 





2.18 





-0.015 


0.775 


1.08 


\i (a.u.) 


12853 




13614 









Table 2. Coefficients C;, in Eq. nSh , for N 2 and NO resonant states. 





Na(X 2 Ii g ) 


NO-( 3 E") 


NO"( 1 A) 


NO"( 1 E+) 


d (eV - 5 "'*) 


0.05 


0.81 


0.60 


0.50 



bottom, has been taken from Ref. |24j . while a is set in order to reproduce the first 
few vibrational levels of Ref. [25]. For NO and NO - the Morse parameters were taken 
from the Refs. [26] [27], [28] . All the constants entering in Eq. ^ for the electronic states 
involved in the processes ([!]) and ^ are given in Table [TJ along with the reduced mass, 
fi, of N 2 and NO molecules. The potential curves, obtained in this manner, are shown 
in Fig.0(a)-(b). 

The width function T(R) in Eq. ^ has been expressed as [T3"] : 

T(R) =Yl Ci \- V ~W ~ V o(R)] kH H (V- - V ) , (8) 

i 

where H is the Heaviside step-function and U is the angular momentum of the lowest 
contributing partial wave associated with the incident electron. The sum in Eq. ^ runs 
over all the electronic states of the molecular negative ion involved in the process, which 
implies one term for N 2 and three for NO - . For the corresponding value of the angular 
quantum number we have l\ — 2 for nitrogen, due to the (i-wave nature of the lowest 
partial wave of the incident electron [15], and l\ — I2 — h — 1 for NO, corresponding to 
a p-wave nature for the three resonant states [21] ■ In order to reproduce the position and 
width of the peaks in the experimental cross section, the constants Ci in Eq. (J8|, for T, 
were considered as phenomenological external parameters of the model and empirically 
adjusted. Their values are given in Table [2] for both the processes Q and @. The 
shift operator A(i?), in Eq. (|5]), has been included in the resonant potential [21]. The 
plots in Fig. [2] show the resonance widths for Ng and NO~, obtained using Eq. (J8|. 
In the case of NO~ the calculated widths are compared with the R-matrix results, 
where the resonance position and width are calculated by fitting the eigenphases sum in 
fixed-nuclei approximation with a Breit-Wigner profile [29]. To validate the numerical 
results, we changed the coefficients in Table [2] by ±10% and noted that the effect on the 
cross sections due to this change is very negligible, in particular for N 2 case. 
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The discrete-continuum coupling matrix element Vk(R) of Eq. ^ is written for 
nitrogen as [21] : 

where k(R) is given by 

2 771 

k 2 (R) = -^[V-(R)-V (R)} . (10) 

For the case of NO we have modulated Vk(R), following Ref. [IS] , with the so-called 
penetration factor, an ad hoc function introduced to force the correct dependence of the 
cross section on the energy near the excitation threshold, namely, 

f(k,R) = l k/k{R) itk<k W (11) 
I 1 otherwise . 

The coupling matrix element for NO is thus written as: 

w-^hm- (12) 



"ilifii — z , m 1X < (13) 



Eq. Q can be related to differential cross section by the approximate expression: 

1 daifip,) 

where the angular factor gi., depending on the scattering angle 9, is explicitly defined in 
Ref. [21] • In the following, for purposes of comparison, we will make use of this equation 
only to extract the integral cross section o"i/(ej) from the measured differential data. 

In the present paper we have computed the cross sections for the transitions 
(v i, Ji) — > (vj, Jf) involving all the possible vibrational levels of N 2 and NO molecules 
and for all the diagonal rotational transitions Ji = Jf = J with < J < 150. 



3. Theoretical and experimental comparisons 

In order to ascertain the accuracy of LCP calculations we compare our resonant 
vibrational excitation (RVE) cross section results for the N 2 and NO molecules with 
the theoretical and experimental cross sections and rate coefficients presently available 
in literature. 



3.1. N 2 

The cross sections for process ([T]) have been calculated by Eq. Q where the degeneracy 
factors for N 2 (X X S+) and N^pIIg) have been set to g — 1 and g r = 4 respectively. 

Fig. [3] compares the cross sections for the rotationless Vi = — > Vf — 1,5,10 
transitions with the experimental data of Allan |9J. The comparison shows a very 
good agreement among the two sets of data for the first two transitions, while some 
discrepancy, around a factor of 2.8 at main peak, is observed in Fig. [3] (b) for the 
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0—7-10 excitation. In this last case, however, the absolute cross section values are 
relatively quite small, which can lead to significant uncertainties in experimental values. 

Fig. [4] compares the theoretical RVE results with unpublished experimental 
measurements of Wong as reported in Ref. [13], for the rotationless vibrational 
transitions — > 1,2. Our calculations are able to reproduce the experimental values of 
the number of peaks, their positions and peak values of cross sections quite accurately. 

Fig. [5] (a)-(b) shows a comparison of the present cross sections with those obtained 
with a non-adiabatic R-matrix approach [H] for the rotationless vibrational transitions 
— > 3 and — > 4 respectively. The relative position of the peaks in the two calculations 
is shifted for both the transitions. This discrepancy can be ascribed to the different 
input quantities used in the two calculations. The RVE cross sections, defined earlier, 
are, in fact, quite sensitive to variations of the potential curves and widths [30J, 
and, as explained in Sec. [2j we have adjusted our input parameters to reproduce the 
experimental data of Allan [9]. 



3.2. NO 

For the RVE of NO molecule, we have calculated the cross sections separately for each 
of the three processes 2 I1 — t- 1 S + , 2 E\ — t- 3 S _ and 2 I1 — t- 1 A. The total cross section has 
been then obtained as a weighted sum of all the contributions, namely [IS], 

fe) = I [< n 4 E+ ^ + ^ + 2 < n 4 A ( e *)] , (14) 

where the coefficients are the statistical weights of each state [18]. Interference terms 
are neglected. 

Fig. [6]shows a comparison of our RVE cross sections for the — > 0, 1, 2, 3 vibrational 
transitions with the differential cross sections measured by Allan [10J at an angle of 
135° and multiplied, according to Ref. [IS] , by a factor of 47r to get, approximately, the 
integrated cross sections. The agreement is quite satisfactory considering the fact that 
we used these experimental results to optimize the input parameters for the adiabatic 
potentials and resonance widths of the Table [T] and |2j No nonresonant background 
contribution has been included in the calculations. In fact, in Fig. [6] (a) the discrepancy 
between our calculated (0 — > 0) elastic cross sections and the measurements of Ref. [TO] 
is precisely due to the contribution of nonresonant background. 

Fig. [7] shows a comparison of our RVE cross sections with the non local calculations 
of Trevisan et al. [18] for the — > 0, 1, 2, 3 vibrational transitions. The agreement here 
also is very satisfactory for all the transitions except for the elastic process, where the 
non local cross sections are significantly higher than the local results. A discrepancy by a 
factor of 6 is in fact observed at the second peak with respect to the present calculations 
which, however, have been calibrated using the experimental data of Allan. 
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3.3. Rate coefficients 

The rate coefficients for both the processes ([I]) and ^ were calculated using Eq. ^ by 
assuming a Maxwellian electron energy distribution function. A comparison for N2 + e 
system, with the well known rate coefficients of Huo |16j . is shown in Fig. [8] for the 
Vi — — > Vf — 1, 2, 3, 4, 5 transitions up to an electronic temperature of T e = 5 eV. The 
representative curves almost overlap, so the two calculations give, in fact, almost the 
same results for the transitions shown in the figure. 

To the best of our knowledge, no rate coefficient data for the electron collision 
processes involving the NO molecule, are available in literature. 



4. Results 



As mentioned earlier, we calculated the cross sections and rate coefficients as a function 
of the electron energy and electron temperature, respectively, for all possible electron- 
impact-induced rovibrational transitions occurring according to the processes ([TJ and 
(j2|. Since the database produced is prohibitively large to be completely shown in the 
present paper, in this Section we limit ourselves to the discussion of some particular 
cases. The whole database, however, will be provided in electronic form elsewhere (see 
Section |5]). 

RVE cross sections for the inelastic rotationless transitions v j = — > Vf = 1 — 5, 10 
for N2 are shown in Figs. [3](5] A general reduction of the magnitude of the main peaks 
is observed as the final level is increased. The same trend is exhibited by the curves of 
Fig. [9] (a)-(b) where the elastic and V{ = 20 — >■ Vf > cross sections, for selected levels 
v f of N 2 , are shown. The trend of the cross sections is reflected in the behavior of the 
rate coefficients, shown in Figs. [9] (c)-(d) as a function of the electron temperature, for 
the same excitation processes as those shown in panels (a)-(b). 

Analogous comments hold for Figs. 10 where the elastic and inelastic cross sections 
as well as rate coefficients for process ([2]), involving the NO molecule, are shown for the 
rotationless vibrational transitions indicated in the figures. 

In view of the plasma modeling applications, we have also investigated the effect of 
the rotational excitation of molecule on cross sections and rates. Fig. 11 (a)-(b) shows 
the cross sections for the two vibrational transitions — > and — > 10, respectively, 
occurring through process ([I]) where now the nitrogen molecule is assumed to be initially 
as well as finally in different rotational levels (Jj = Jf = J). We have in particular 
considered the three cases where J = 0, 100 and 150. No significant variation is 
observed in the order of magnitude of the cross sections in going from J = to J = 
150. An enhancement of the main peaks, by about a factor of 2, is seen in fact for the 



vibrationally elastic process in Fig. 11 (a), while a decreasing of the peak values occurs 



for the vibrationally inelastic transition in Fig. 11 (b) as J is increased. A shift toward 



smaller energy is also observed, as expected, in both the excitations, coming from the 
lowering of the threshold caused by a reduction of the transition energy. Fig. 12 (a)-(b) 
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shows the same cross sections occurring through process ([2]) for NO molecule. However, 
in this case the main peaks are reduced for vibrationally elastic process and enhanced 
for vibrationally inelastic process, with respect to J = 0, as the value of J is increased. 



NO, respectively, for different J values and for the elastic — > and inelastic — > 10 
vibrational transitions. Inspection of the figures for the elastic case shows that the 
rates are more sensitive to the rotational contribution for temperatures above 1 eV for 
nitrogen and above 0.2 eV for nitric oxide. For T e < 1 eV practically no rotational effects 
is observed in N 2 , while for NO the rates become insensitive to the initial rotational 
excitation of the molecule for high J and for T e < 0.2 eV. For the inelastic case the 
situation is less clear, but no large variations, as J changes, are observed in the whole 
range of temperatures investigated. 

5. Summary 

We presented in this paper the vibrational excitation cross sections and rate coefficients 
for processes occurring through resonant electron-impact on N 2 and NO molecules, 
for modeling purposes. Complete sets of data have been produced by using the local 
complex potential model for the description of the nuclear dynamics, and by representing 
the required input electronic potential energies and widths by parametric standard 
expressions, appropriately adjusted to reproduce the experimental measurements. We 
also illustrated the calculation by commenting on some particular selected Vi — > Vf 
transition, and showing the effect of the initial rotational excitation of the target 
molecule on RVE cross sections and rates. 

The entire cross section and rate coefficient database, including all the — > 
Vf transitions, and for J ranging from through 150 is available at the website: 
http:/ /users. ba.cnr.it/imip/cscpal38/phys4entry/database. html 
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Figure 1. Morse potentials for: (a) N2 (solid line) and N^~ (dashed line); (b) NO 
(solid curve) and the three NO - states: 3 E~ (short dashed line), *A (dot-dashed line) 
and 1 S + (long dashed line). 
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Figure 3. Electron-N2 theoretical-experimental cross section comparison for the 
rotationless Vi = — > Vf = 1,5, 10 transitions as indicated in the panels. Full lines: 
present calculations; dashed lines: experimental results of Allan [9]. 




1 2 3 4 51 2 3 4 



Electron energy (eV) Electron energy (eV) 



Figure 4. Electron-N2 theoretical-experimental cross section comparison for the 
rotationless Vi = — > Vf = 1,2 transitions as indicated in the panels. Full lines: 
present calculations; dashed lines: unpublished experimental results of Wong as cited 
in Ref. [13]. 
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Figure 5. Electron-N2 comparison of the present RVE calculations (full line) with the 
R-matrix results (dashed line) [14] for the rotationless transitions shown in the panels. 




Figure 6. Electron-NO comparison of the present RVE calculations (full lines) with 
the measurements of ref. [TU] (dashed lines) multiplied by a factor of 4ir (see text) for 
the vibrational transitions — > 0, 1, 2, 3. 
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Figure 7. Electron-NO comparison of the present RVE calculations (full lines) with 
those of Trevisan et al. [IS (dashed lines) for the vibrational transitions —> 0, 1, 2, 3. 




Figure 8. Comparison of the present electron-N2 RVE rate coefficients (full lines) 
with those of Huo et al. [16] (dashed lines) as a function of the electron temperature. 
The curves refer to the transitions starting from Vi = and ending on the final levels 
Vf = 1, 2, 3,4, 5 as indicated in the figure. 
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Figure 9. N2 vibrational excitation cross sections and corresponding rate coefficients 
for (a, c) Vi = Vf and (b, d) Vi — 20 — > Vf > Vi. 
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Figure 11. Cross sections and corresponding rate coefficients for (a, c) — > elastic 
and (b, d) — > 10 inelastic vibrational transitions involving the N2 molecule initially 
in different rotational levels: full line J = 0, dashed line J = 100, dotted line J = 150. 
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Figure 12. Same as Fig. [IT] for c-NO process. 
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